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Abstract 

We present numerical data of the height-height correlation function and of the 
avalanche size distribution for the three dimensional Toom interface. The height- 
height correlation function behaves samely as the interfacial fluctuation width, 
which diverges logarithmically with space and time for both unbiased and biased 
cases. The avalanche size defined by the number of changing sites caused by 
a single noise process, exhibits an exponentially decaying distribution, which is 
in contrast to power-law distributions appearing in typical self-organized critical 
phenomena. We also generalize the Toom model into arbitrary dimensions. 
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The Toom model [1] is a dynamical model of Ising spins in which the condition 
of detailed balance is not satisfied and hence whose phases are not described 
by equilibrium Gibbs ensembles. Recently Derrida, Lebowitz, Speer and Spohn 
(DLSS) [2] studied physical properties of interfaces formed in the two-dimensional 
Toom model. In low- noise limit, this model leads to a (1 + 1) dimensional solid- 
on-solid-type (SOS) model, which is in turn much simpler for understanding 
generic nature of dynamics. In the SOS model, the dynamics of spin-flips may be 
regarded as a deposition-evaporation process of particles. Due to the nature of 
the Toom dynamics, the deposition-evaporation process occurs in an avalanche 
fashion with preferred direction. DLSS found, among others, that the continuum 
stochastic equation describing such SOS model is the well-known Kardar-Parisi- 
Zhang (KPZ) equation [3]. In reference [4], present authors proposed a natural 
generalization of the model in three dimension on the body centered cubic lattice 
and found that the continuum equation is the anisotropic KPZ (AKPZ) equation 
given by 

d t h = u^h + u ± d ± 2 h + \\\\{d\\h) 2 + h±(d ± h) 2 + 77 (1) 

with opposite signs of the coefficients Ay and Aj_. In Eq. (1), h is the height of the 
fluctuating surface with respect to a reference (co-moving) plane, d\\ (d±) stands 
for spatial derivative along the direction parallel (perpendicular) to the avalanche 
direction, and 77 is the white noise. AKPZ equation with opposite signs of d» and 
<9_i_ is known to renormalize toward the weak coupling limit [5] and consequently, 
square of the width, w 2 , of the Toom interface shows logarithmic dependence 
both in space and time. This is in contrast to the model studied by Barabasi, 
Araujo, and Stanley [6] which again is described by Eq. (1) but belongs to the 
strong coupling regime of the AKPZ universality. 
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In this work, we present numerical data of the height-height correlation func- 
tion and of the avalanche size distribution for the three dimensional Toom inter- 
face. The height-height correlation function shows the same logarithmic depen- 
dences both on space and time as in w 2 . The avalanche size distribution is found 
to be exponential. This is in contrast to power law behaviors appearing in typical 
self-organized critical phenomena and indirectly substantiate the validity of the 
collective variable approximation used by DLSS. 

The d-dimensional Toom model we introduced in [4] and generalized here con- 
sists of Ising spins (a(xi,X2, • • • , Xd) = ±1) on rf-dimensional body-centered-cubic 
lattice. The spin coordinates Xi takes the values in Z for the spins on one sublat- 
tice and Z+| for those on the other sublattice. At each time step, a randomly 
selected spin is updated according to the local rule that it becomes, at the next 
time step, equal to the majority of itself and of its 2 d ~ 1 neighbors relatively situ- 
ated at (— 1,±|,±|,---,±|) with probability 1—p — q, to +1 with probability p, 
and to —1 with probability q. The referencing neighbors are shown in Fig. 1 for 
d — 3. Unbiased (biased) dynamics results when p — q {p ^ q). To produce sta- 
tionary interface at zero noise, the boundary spins on the three surfaces defined 
by x\ = 0, X2 = and X2 = L, respectively, need to be fixed to the value +1 (— 1) 
if x 2 > (<) L/2 where L is the system size. Periodic boundary conditions are 
imposed to other surfaces. 

In the low- noise limit, a spin flipped due to noise returns immediately to its 
original state by the majority rule, if the spin is situated away from the inter- 
face. Accordingly, we may assume that spin flips occur only at the interface. We 
are thus led to study an effective SOS-type model as the generalization of the 
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(1 + 1) -dimensional stairlike model studied by DLSS. In the SOS model picture, 
the Toom dynamics is mapped to a particle dynamics on a d s = d — 1 dimensional 
substrate in the form of deposition-evaporation with avalanche. The (d s + 1) di- 
mensional SOS model is defined on d s dimensional body-centered-cubic lattice. 
For d s = 2, it is the checkerboard lattice, square lattice rotated by 45 degrees. 
To each lattice point, a relative height corresponding to the X2 coordinate of the 
original Toom interface is assigned. The height should satisfy the SOS condition, 
that nearest neighbor heights should differ only by ±1. Initially we begin with a 
flat surface characterized by height on one sublattice and height 1 on the other. 
At each time step, we select a random site and start evaporation (deposition) 
process with probability p (probability 1 — p) . In the evaporation (deposition) 
process, the height of the selected site is decreased (increased) by 2 if the SOS 
condition is satisfied with respect to the 2 da ~ 1 spins to the left. By left spins, we 
mean those nearest neighbor spins whose x\ coordinate is less by \. If the SOS 
condition to the left is not satisfied there is no change. Next the avalanche pro- 
cess proceeds to the right (+xi direction) until all sites satisfy the SOS condition. 

When d s — 1, this avalanche dynamics reduces to the spin exchange dynamics 
of DLSS. The unbiased (biased) dynamics corresponds to the case p — \ {p ^ \). 
If the avalanche process is not allowed, so that depositions and evaporations oc- 
cur only on local valleys and mountains respectively, then our model would be 
equivalent to the deposition-evaporation model proposed by Forrest and Tang [7] , 
a generalization of the Plischke-Racz-Liu model [8] into higher dimensions. More 
detailed description of the SOS model for d s = 2 is found in the original paper 
[4]. Although the model is defined for general d — d s + l, we confine our attention 
to d = 3 from now on. 
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The continuum equation for the (2 + l)-dimensional model is described by 
Eq. (1), because it selects out a preferred direction, and because the cubic non- 
linear term derived by DLSS was proved to be marginally irrelevant even in (1 + 1) 
dimensions [9]. For the unbiased case, both of the nonlinear terms in Eq. (1) dis- 
appear due to the symmetry of deposition and evaporation. So the equation 
reduces to the Edwards- Wilkinson (EW) equation [10] implying that the square 
of the surface width diverges logarithmically with space and time. For the biased 
case, average height grows with increasing time, so that Ay and Aj_ are nonzero. 
It was shown [4] that Ay > and Aj_ < by applying the tilt argument [11]. 
Consequently, the model belongs to the weak-coupling regime of the AKPZ uni- 
versality. Therefore, both the square of the surface width and the height-height 
correlation function are logarithmic for both unbiased and biased cases. 

To check correspondences between the low noise Toom model and the SOS 
model, we performed numerical simulations for both models. Results for both 
cases run on small sizes are in complete agreement with each other in low-noise 
limit. Accordingly, we performed simulations intensively for larger systems using 
the SOS model. The simulations are done in the range of system size L = 20 ~ 
140 for both unbiased (p = 0.5) and biased (p ^ 0.5) cases. We measured 
the surface width w 2 and the height-height correlation functions C\\(r\\,t) and 
C±(r±,t), defined by 

C\\(r h t) = {[h(r h t)-h(0,0)] 2 ), and C ± (r ± , t) = ([h(r ± , t) - h(0, 0)] 2 ), (2) 

respectively. Here, ry (r±) denotes the spatial separation along the direction 
parallel (perpendicular) to the avalanche direction. From the theoretical con- 
siderations, we expect the height-height correlation functions diverges as ~ lnt 
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before saturation, and ~ lnr after saturation. Typical data of C\\(r) and C±(r) 
after saturation are shown in Fig. 2 for the case of unbiased dynamics (p = 0.5) 
and in Fig. 3 for the case of biased one (p = 0.3). Each curve is for system sizes 
L = 40, 60, 80, 100, and 120, respectively, from bottom to top and is averaged 
over 300 configurations. It confirms clearly the logarithmic dependence on r. 
Deviation from straight line for large r is the finite size effect. 

Next, we examined the avalanche size distribution n(s). The avalanche size s 
is defined as the number of successive spin flips by a single noise process. n(s) 
was measured in two different manners. In the first case, it is measured in the 
critical state (after saturation), while in the second case, it is measured during 
the whole time steps. In both cases, the distribution function n(s) is found to 
be exponential, n(s) ~ exp(— s/s*) as shown in Figs. 4 and 5. The characteristic 
size s* determined from the slopes of the inset figures is found to be independent 
of p and of the way n(s) is measured. It takes the value s* = 1.11 ± 0.01. 

In conclusion, we have generalized the Toom model into d dimensions and de- 
fined its associated SOS-type model. Also we have presented numerical data for 
the height-height correlation functions in perpendicular and parallel directions, 
respectively, for unbiased and biased cases for d — 3. We have found that for the 
unbiased case, the interface is described by the EW equation, and for the biased 
case, it is described by the AKPZ equation with the opposite signs of Ay and 
Aj_. Consequently, the height-height correlation functions diverge logarithmically 
with space and time. The avalanche size distribution has also been examined, 
which exhibits an exponential distribution for unbiased and biased cases instead 
of a power-law distribution. 
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Figure Captions 

Fig. 1 The three dimensional Toom rule on bcc lattice used in this work. The 
black circled spin in the center is updated with the majority rule of itself 
and four nearest neighbor spins (the black circled) with probability 1— p — q, 
and becomes equal to +1 (—1) with probability p (q). In the low-noise limit, 
p,q->0. 

Fig. 2 The height-height correlation functions C\\(r) (dotted line) and C±(r) 
(solid line) versus lnr for unbiased case, after saturation. Each curve is for 
system sizes L = 40, 60, 80, 100 and 120, respectively, from bottom to top, 
and is averaged over 300 configurations. 

Fig. 3 Same as in Fig. 2 for biased case, p = 0.3. 

Fig. 4 The avalanche size distribution n(s) versus s for unbiased case measured 
during whole time steps (a) and after saturation (b). The data are obtained 
from different system size, but they are collapsed into each other, implying 
that n(s) is independent of system size. Insets show lnra(s) versus s. 

Fig. 5 Same as in Fig. 4 for biased case, p = 0.3. 
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